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Abstract 



Quantum-dot Cellular Automata ( QCA ) provides a basis for classical computation without tran- 
sistors. Many simulations of QCA rely upon the so-called Intercellular Hartree Approximation 
(ICHA), which neglects the possibility of entanglement between cells. Here, we present computa- 
tional results that treat small groups of QCA cells with a Hamiltonian analogous to a quantum 
mechanical Ising-like spin chain in a transverse field, including the effects of intercellular entan- 
glement. When energy relaxation is included in the model, we find that intercellular entanglement 
changes the qualitative behaviour of the system, and new features appear. In clocked QCA, isolated 
groups of active cells experience oscillations in their polarization states as information propagates. 
Additionally, energy relaxation tends to bring groups of cells to an unpolarized ground state. This 
contrasts with the results of previous simulations which employed the ICHA. The ICHA is a valid 
approximation in the limit of very low tunneling rates, which can be realized in lithographically 
defined quantum- dots . However, in molecular and atomic implementations of QCA, entangle- 
ment will play a greater role. The degree to which entanglement poses a problem for memory and 
clocking depends upon the interaction of the system with its environment, as well as the system 's 
internal dynamics. 

1 Introduction 

Quantum-dot cellular automata (QCA) provides a basis for binary computation that is fun- 
damentally different from today's transistor-based technology [TJ[2] . With QCA, information is 
stored and transferred by cells of quantum dots, such as those depicted in Figure [1} A cell consists 
of four quantum dots arranged at the corners of a square. Each cell contains two mobile electrons, 
which can tunnel between dots. The electrons experience mutual repulsion and at sufficiently 
low temperature will take either of the two possible antipodal configurations (or a superposition 



thereof). Because of the symmetry of the cell, neither configuration is preferred in the absence of 
a perturbation, and the cell, with its two electrons, can be thought of as a bistable switch: a small 
push toward one or the other diagonal will cause the electrons to "switch" almost completely. 
The interaction between cells tends to align the electronic configuration of one cell with that of its 
neighbours. A line of cells thus acts as a binary wire. Different geometrical arrangements of cells 
correspond to logic gates, and together enable the design of a universal computer capable of very 
low power operation |3j|4] . Computing speed would be enormously improved by the incorporation 
of clocking (5j|6), whereby clock zones are used to sequentially send "bit packets" through a QCA 
array to allow for pipelining. A clock zone consists of a group of cells whose parameters can 
be independently tuned to either allow or disallow the group of cells to become polarized (i.e., 
to take on one of the antipodal configurations). A bit packet consists of a group of interacting 
cells in "active" clock zones. As envisioned theoretically, many bit packets could be processed 
simultaneously in a single QCA circuit by modulating clock zone parameters appropriately. 

QCA cells were demonstrated experimentally in metal-island quantum dots as early as 1997 [7]. 
These metal-island QCA cells have been used as prototypical devices to demonstrate transmission 



of information, logic gates |8l 9], and clocking 10,11 . In this implementation, clocking is realized 



by modulating the null-dot potential between dots. However, the relatively large size of the 
quantum dots means that energy levels are closely spaced, so metal-island devices must be kept 



at temperatures below ~ 5K 12 13 in order for quantum effects to be observable. As the 
dimensions of a QCA cell are reduced, the operating temperature increases, and at the molecular 
scale room temperature operation becomes possible. Further advantages of miniaturization include 
fast switching times and increased device density. For these reasons, molecular scale QCA has 



held great promise. Suitable candidates have been synthesized 14, 15 , and a QCA cell made of 
atomic quantum dots on silicon has also been realized (16). For QCA cells of this size it is probably 
neither feasible nor desirable to have addressable control over the parameters of individual cells. 
Instead, large groups of cells could be addressed by external fields (5 , 17, 18 . The variation in 
these external fields can produce clock zones. 

A great deal of theoretical and modelling work has been done on the topic of QCA [2 19- 



23 . This research has aimed to capture the qualitative and quantitative characteristics of QCA 
cells and of arrays of cells. Because of the difficulties inherent in solving the complete quantum 
mechanical problem, a number of simplifying assumptions are typically made. These include 



a reduction of the Hilbert space to two states per cell 20 , treatment of intercell interactions 
via a mean- field approach [2, 19 , and finally an assumption of exponential energy relaxation [3j 
21 . Although much insight has already been gained on the nature and importance of quantum 



mechanical calculations that go beyond these approximations, including full and partial quantum 



correlations 20,21 , some important features of clocked QCA systems remain unreported. We 
have found that the inclusion of intercell correlations can significantly change the steady state of 
the system not only quantitatively, but qualitatively, even in the case of very simple systems, such 
as an unbiased line of QCA cells. Full quantum mechanical simulations, with the approximation 
of exponential relaxation to a thermal steady-state, predict an exponential loss of information 
even as bit packets propagate, as well as coherent oscillations whose period strongly depends 
on the size of the bit packet. This contrasts with mean-field simulations, which show indefinite 
propagation of information. The results presented here will have implications for molecular-scale 
QCA device design, and will highlight the need for implementation-specific theoretical treatments 
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of the interaction of a QCA system with its environment. 

In Section [2] we review the basic theory of QCA as well as some of the most common approx- 
imations used in QCA simulations. In Section [3] we consider the full Hamiltonian for a QCA 
line, and the characteristics of its solutions. In Section [4] we present the results of fully quantum 
mechanical simulations. Section [5] discusses the main results of the paper, their scope, and their 
implications for QCA design, and finally we offer a conclusion in Section [6] 

2 Simulation of QCA Systems 



The dynamic behaviour of QCA was first explored in 20 , where Tougaw et al. examined the 
time evolution of QCA cells by considering a basis set consisting of all sixteen possible states of 
a four-dot QCA cell populated by two electrons of opposite spin. In this sixteen-state approach, 



the authors construct a Hubbard-type Hamiltonian, given as 20 
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where the operator dj jCr (m) (a] CT (m)) annihilates (creates) an electron on the i th site of cell m 
with spin a, the operator nj jCr (m) = aJ CT (m)dj i(T (m) is the number operator for an electron on 
the i th site of cell m with spin a, and Vq = ql/(4:7ie) is a constant where q e is the charge of 
the electron and e the electrical permittivity of the medium. The first term in Eq. [T] represents 
the on-site energy of a dot. The second term describes the electron tunnelling, where ty is a 
hopping constant (with units of energy) between neighbouring sites i and j, determined from the 
structure of the potential barriers between the dots in the cell. The third term in Eq. [T] accounts 
for the energetic cost, Eq, of putting two electrons of opposite spin at the same site, and the final 
two terms are related to the Coulombic interactions between electrons in the same cell and in 
neighbouring cells, respectively. The polarization of each cell can then be found by evaluating 

_ (p? + pf) - (pi - p?) 

PT + P2+PT + PT ' 

where p™ is the expectation value of the number operator on the i th site of cell m; i. e., p™ = (hi{m)). 
Sites within a cell are labeled clockwise starting from the top right, as shown in Figure [TJ While 
this Hamiltonian considers the complete many-body configuration space, including correlation ef- 
fects within and between cells, it becomes computationally intractable when used to model large 
systems. For example, a three cell system requires 16 3 = 4096 basis vectors. The exponential 
growth of the basis set makes it computationally prohibitive to model any circuit larger than just 
a few cells. We will now discuss three ubiquitous approximations used for studying QCA circuits 
and systems. 
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2.1 Two-State Approximation 



As shown in 20 , the ground state of a single cell within the full sixteen-dimensional Hilbert 
space remains almost completely contained within the two-dimensional subspace of the polarized 
basis vectors. For sufficiently low temperatures, we can therefore expect the state of a cell or a 
line to be well described by the two-state approximation. We refer to this reduced basis as the 
polarization basis, and denote the two states as |0) and |1), as shown in Figure [TJ The Hamiltonian 
in the polarization basis for a system of N interacting QCA cells, under the influence of driver 
cells, is then described by a 2 N x 2 N Ising-like Hamiltonian: 



JV -. N 1 N 

H = - J2 - 2 E + 2 E E tffPrtM ( 3 ) 

i=l i<j D 8=1 

where 7* is an effective tunneling energy, related to the hopping energy, Uj, in equation [TJ E]f 
is the so-called kink energy between cells i and j, and accounts for the energetic cost of two cells 
having opposite polarization. Pjj labels the polarization of the driver cell labelled D, and E]f* is 
the kink energy between cell i and driver D. The driver cells provide a mechanism for input into 
the QCA circuit and have a polarization that can range from —1 to +1. The Pauli operators for 
the i th cell, & a (i); a — x,y, z, represent the tensor product of iV 2 x 2 identity operators, with the 
i th identity operator replaced by the one of the Pauli matrices, 

1 \ A / i \ A / -1 



1 1 j' ay \ -i J ' az v 1 

For example, <r y (2) ^l®^®!®...®!. The polarization of a cell, i, can now be defined as 
Pi = —(a z (i)). The first term in Eq. [3] accounts for the kinetic energy of electrons, and tends to 
bring the cells to a superposition of polarization states. The second and third terms account for 
the energy cost of having a cell misaligned with its neighbour, or with a driver cell respectively. 
The Hamiltonian in Eq. [3] is very similar to the quantum Ising model with a transverse field, 



which has been extensively studied 24 , 25 



In Section [3] we will use the two-state approximation to study systems of interacting QCA cells. 
This procedure, like the full sixteen-state one, naturally includes the effects of inter-cell entangle- 
ment, or quantum correlations. The price we pay for proper inclusion of quantum correlations is 
an exponential growth in the basis set with the number of interacting cells, N. Even with the 
two-state approximation, this limits its application to systems containing only a small number of 
interacting cells. Further approximations are used to solve the problem of scaling. 

2.2 Intercellular Hartree Approximation 

One way to eliminate the problem of exponential scaling is to ignore inter-cellular entanglement 
effects altogether and solve the Schrodinger equation for each individual cell separately. This 
method is known as the intercellular Hartree approximation (ICHA) (2|[l9] . In this Hartree-type 
treatment, cells are coupled to one another via expectation values (polarizations) rather than 
operators. The Hamiltonian (in the polarization basis) for a single cell i, is then simply 
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Figure 1. Schematic representation of the two diagonal states which form the basis in the two-state 
approximation. 



Hi = -Ji(T x + 



(4) 



where Pj is the polarization of cell j. The polarization of cell i is found by evaluating, Pj = — (er z ). 
Because the solutions of one cell's Hamiltonian define parameters that enter its neighbours' Hamil- 
tonians, the system of Schrodinger equations must be solved iteratively to obtain self-consistency. 
In calculating the state at a particular time, the initial guess is typically taken to be the state of 
the system at the previous time step. 

The primary benefit of this approximation is that it only requires the diagonalization of iV 2 x 2 
Hamiltonians, which scales linearly with the number of cells in the system. Furthermore, if 
a circuit evolves adiabatically and cells are assumed to remain at the ground state, then this 
problem simplifies further by recognizing that the polarization of any cell, i, can be evaluated 



analytically using 26 



1 E4 J >; 



27 



P 



\ 



(5) 



1 + 



Eq. [5] produces the well-known nonlinear cell-to-cell response function shown in Figure [2] While 
the ICHA is generally capable of arriving at the correct ground state of an array of QCA cells in a 
single clocking zone with a fixed driver, it can falsely predict a latching mechanism within a group 
of cells that allows them to retain (and even obtain) a polarization in the absence of a perturbing 
cell. This goes beyond previously noted inaccuracies of the ICHA, such as in dynamics and finite 
temperature behaviour |20|, where the many-cell excited states are needed to get quantitatively 
correct results. 

To illustrate the effects of the ICHA on the calculated ground state of QCA arrays, consider the 
driven two-cell wire shown in Figure |3](a). Two simulations were conducted on the wire; the first 
using the ICHA, and a second using the more complete quantum mechanical treatment discussed 
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Figure 2. Nonlinear cell-to-cell response function for 7 -C E^ and ksT <C (optimal conditions 
for QCA). The output cell is almost completely polarized for even a small driver polarization. 



in Section [2T For each simulation, the driver cell polarization was swept between, — 1 < P < +1, 
and repeated for cycles with different values of E%./j. The polarization of the first cell in each of 
these simulations is plotted in Figures [3|b)-(c). 

Let us first consider the results from the ICHA simulation. When using the ICHA, one has 
freedom in choosing an initial guess for the polarizations. The most common method in dynamic 
simulations is to use as an initial guess at each time step the solution from the previous time 
step. With this method, Cell l's response to the driver cell follows the hysteresis curve shown in 
Figure |3](b). There is a retained polarization even as the driver polarization goes to zero, which 
is indicative of a memory effect. This type of hysteresis curve is commonly seen in ferromagnets. 
When a magnetic field is applied to a ferromagnet, electron spins in the ferromagnetic material 
align themselves with it. For temperatures below the Curie temperature, this alignment is retained 
when the field is removed, and the material remains (at least partially) magnetized until a magnetic 
field in the opposite direction is applied. The ICHA predicts an analogous behaviour for QCA 
cells, with the role of the external magnetic field, in this case, played by the driver cell. When the 
driver cell is turned on, both the driven cells align themselves with the driver polarization. As the 
driver cell polarization is removed, the coupling between the two driven cells (through expectation 
values) allows them to retain part of their polarization even as the driver cell polarization reaches 
zero. Only after a sufficiently strong driver polarization in the opposite direction do both cells 
switch polarizations. The strength of the driver polarization required to switch the driven cells 
depends solely on the ratio j/E^', the lower the tunnelling rate, the larger the residual polarization 
of the driven cells as the driver cell's polarization is removed. 

The other method that can be used with the ICHA is to sample the space of polarizations in 
search of the self-consistent solution with the lowest energy. For certain driver cell polarizations, 
there are two self-consistent solutions. In Figure ^b), the lowest energy solution is shown in bold. 
If the lowest energy solution is always used, the ICHA predicts a discontinuity in the response of 
the driven cells, at zero driver cell polarization. As the driver's polarization crosses zero, the cells 
respond by abruptly "snapping" to the other polarization state. The ICHA is not typically used 
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Figure 3. Parametric simulations of a driven two-cell wire for different tunnelling rates, (a) The 
two-cell wire being simulated, (b) Polarization of Cell 1 as a function of the driver cell polarization, 
calculated using the ICHA. The ICHA predicts a hysteretic response indicating a memory effect. 
Where two solutions exist, the one with lowest energy is shown in bold, (c) Polarization of Cell 1 as 
a function of the driver cell polarization, calculated using the more complete quantum mechanical 
treatment described by equation [3] This more complete simulation shows no hysteresis. We have 
taken k B T « E k /A, which corresponds to E k w 100 meV for room temperature. 



in this way, however. 

Simulations conducted using a more complete quantum mechanical treatment show no such 
hysteresis, as shown in Figure |3]^c). Here, as the driver polarization is removed, Cell 1 also relaxes 
to zero polarization. Thus, there is a "depolarizing" effect that is not predicted when treating 
QCA systems using the ICHA. Also, the response is continuous, but still non-linear and with a 
slope greater than unity at the origin. 

The ICHA represents, in a sense, the minimum inclusion of quantum mechanical effects and 
it will be shown that this assumption can lead to errors when applied to clocked QCA. It is 
possible to make a less drastic simplification of the system dynamics by including some, but not 



all, intercell correlations, for example by including only nearest-neighbour pair correlations 21 22 



Including some of the correlations has been shown to give the correct answer in some cases where 
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the ICHA fails. In Section [4] we will show a case in which the ICHA fails, and the inclusion of 
nearest neighbour correlations does not help. 

2.3 Relaxation Time Approximation 

A final approximation to be discussed here is the relaxation time approximation, which is 
commonly used in simulations of QCA dynamics. In the absence of energy dissipation and other 
decohering effects, a QCA array will evolve coherently according to the Liouville equation for the 
density matrix, p: 

(6) 

which, for a pure state, is exactly equivalent to the Schrodinger equation. The Hamiltonian in 
Eq. [6] may be time-dependent, e.g., when driver cells are switched or when cell parameters are 
changed to implement clocking. 

Over fairly short time scales, quantum mechanical systems often fall to a thermal steady 
state (3]. If the QCA system is weakly-coupled to the environment, and the energy transfer 
between the system and environment is well-described by a Markov process, then at low temper- 
atures, the simplest way to incorporate energy dissipation into a model of QCA dynamics is via 
the relaxation time approximation [3j[2TJ[27j[28]. This is done by adding a dissipation term to 
Eq.[6j 

" " (P(t) - Pss) , (7) 
T 

where r is a phenomenological time constant, and p ss is the steady-state matrix defined as 

e -H(t)/k B T 

Pss = 7 ; =r- (8) 

Tr <^ e-n(t)lk B T I 

Determining the steady state density matrix exactly is tantamount to solving the complete 
Schrodinger equation for the system, and thus comes up against all the difficulties mentioned 
above. It is therefore usually calculated using the ICHA and the two-state approximation. 

Eq. [7] imposes an exponential approach of the density matrix towards its steady-state value, 
with a time constant r. This is one example of a quantum master equation, or an equation of 
motion for the density matrix, which in this case is phenomenological. It would be useful to derive 
a quantum master equation based on the microscopic details of the system, its environment, and 
their interaction. This would show whether or not the form of Eq. [7] is correct, and how the 
value of t relates to microscopic parameters. However, a precise quantum master equation will 
in general be implementation dependent. In this paper, we treat an idealized QCA that is not 
tied to a specific implementation. Furthermore, regardless of the specifics of the quantum master 
equation, the density matrix in Eq. [8] will very often represent the real steady state of atomic 
and molecular systems. We therefore attempt to calculate the correct steady state behaviour, 
acknowledging that the dynamics and the specific value of r will be implementation-dependent. 



dt. 



Pit) 



1 

ih 



H(t),p(t) 



dt 



Pit) 



Hi 



Hit), pit) 



8 



3 Full Quantum Mechanical Calculations 



Having developed an appreciation for the effects of the ICHA and of the relaxation time ap- 
proximation on the calculated dynamics of QCA arrays, we will briefly consider the qualitative 
features of some tractable systems of interacting QCA cells by once again considering the two-state 
Hamiltonian described in equation |3j 

The difficulty in solving Eq. [3] can be reduced by assuming only nearest-neighbour coupling so 
that E]f = Ek5i±\j and E^ D = for cells that are not adjacent to the driver labelled D. This 
does not imply that correlations beyond nearest neighbours are ignored, however. We also restrict 
our attention to linear chains of cells. 
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(b) Analogous Double Well Systems 

Figure 4. (a) Spectra of unbiased lines of interacting QCA cells, ranging in length from one to eight 
cells, with j/E k = 0.1. The spectra are solutions of Eq.[3j with a constant (N — l)E k /2 added 
for ease of interpretation. The inset shows the difference in energy between the two lowest energy 
levels for each line, (b) Conceptually, the array can be viewed as a single two state system with a 
barrier that increases with the number of cells. 



The spectrum for an unbiased line of N (ranging from 1 to 8) QCA cells, within the two-state 
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approximation is shown in Figure |4^a). For a line of N cells, there are 2 N eigenstates. Under the 
conditions required for QCA-based computing, that is 7 << E^, the energy levels come in clusters 
separated roughly by E^. There are always two non-degenerate lowest-energy states. In the case 
of a single cell, these are in fact the only two states, and their separation is exactly 27. They 
correspond to the symmetric (\ip s ) — + aric ^ anti-symmetric (\ip a ) — ^(1^) ~~ |1))) 

combinations of the polarization states (with the former being the ground state). We note that 
the polarization basis vectors, |0) and |1), are not energy eigenstates of the undriven cell. If a 
single quantum measurement of a z is carried out on a QCA cell in its ground state (or in its first 
excited state for that matter), the outcome will yield either —1 or +1, with equal probability, 
i.e., P = 0. In this sense, one can say that the ground state of a single unbiased cell carries no 
information. Only in the limit where 7 — > do the polarization basis vectors become valid energy 
eigenstates. 

For longer lines, the separation between the two lowest-lying energy eigenstates becomes smaller 
and they represent entangled states. Specifically, they are the symmetric and anti-symmetric 
combinations of the state with all cell polarizations aligned along one diagonal, and the state with 
them aligned along the other. That is, \if) s>a ) ~ ^(|000 . . . 0) ± 1 1 1 1 . . . 1)). The equality becomes 
exact in the limit where Ej t / r y — > 00. 

Again, we note that the "aligned" states, 1 000 ... 0) and | 111 ... 1), are not the energy eigen- 
states. The energy eigenstates are in fact superpositions of the aligned states, and the polarization 
of any cell in a line in its ground state (or first excited state) is P = 0. This indicates that the 
eigenstates do not carry information. The inset of Figure |4|a) shows that the splitting between 
these two lowest-lying states decreases exponentially. Specifically, each added cell causes the split- 
ting to decrease by a factor of ~ Ek/j. Only in the limit where this splitting goes to zero (i.e. an 
infinite number of interacting cells), do the aligned states become energy eigenstates [25). With 
the decrease in the splitting of the two lowest energy levels, comes an increase in the time required 
for coherent tunnelling from one polarization state to the other. For low temperatures, we can 
think of the group of cells as a two-state system, analogous to a double well, as depicted in Figure 
P](b). As the length of the line increases, the barrier separating the two aligned states increases. 
It follows that longer lines acquire increased bistability. However, at finite temperature, thermal 
fluctuations may cause excitations and an eventual approach to the unpolarized steady state. 

Figure [5] shows the unitary time evolution of unperturbed lines of one, three, and five cells 
initially in an aligned state in each case. All the cells in the wire oscillate in unison between the 
two polarization states. This simulation represents the limit of infinite relaxation time (i.e., no 
energy dissipation or loss of phase coherence). Since the longer lines exhibit increased bistability, 
the polarization state can be maintained for an arbitrary length of time by increasing the number 
of cells in the line. For single cells and small groups at the atomic scale, coherent oscillations 
will likely be much faster than the measurement time, and therefore will lead to a loss of classical 
information. For larger groups, the loss of information will likely be limited by environment- 
induced phase decoherence and energy relaxation. 

Finally, we emphasize that short unbiased lines, like individual cells, have a unique ground 
state. This contrasts with the commonly made assertion that cells and lines have two degenerate 



states 29 , namely, the aligned states. The ground state is an entangled state, and is therefore not 
accessible to the ICHA. Based on this observation, as well as the above-mentioned shortcomings 
of the ICHA in predicting the behaviour of even very short lines, we are led to a more in-depth 
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investigation of QCA line dynamics, particularly with regard to clock zones which are spatially 
separated from any driver cells. 




Figure 5. Coherent oscillations of the first cell in unperturbed one, three, and five cell lines. The cells 
are all initially in the P = 1 state. In each line, the cells oscillate together, so the polarizations of the 
first cell gives a good representation of the polarization of all its neighbours. The fast oscillations in 
the 3 and 5 cell lines are due to a small component of higher-energy states, which manifest as kinks 
propagating and reflecting through the line. 7 = 10 meV and E k = 100 meV. 



4 Loss of Polarization in Isolated Bit Packets 



As discussed in Section 2.2 the ICHA predicts a latching mechanism within a line of interacting 
cells that allows them to polarize (and retain this polarization) in the absence of a fixed driver 
cell. This phenomenon is found to be an artifact of the ICHA, not necessarily representative of 
actual dynamics. Consider a single, unbiased QCA cell. At low temperatures, the QCA cell will 
relax to its ground state, \ip s ) = 77g(|0) + and thus will have a polarization, P = 0. At higher 
temperatures, interactions with the environment can cause the QCA cell to be in a mixed state, 
and it is best described by a density matrix. If the steady state density matrix of the system 
is the one described in Eq. [8j then it is easy to show that the expectation value of o z is zero, 
and therefore the polarization is also zero. For longer lines, the same reasoning applies, the only 
differences being that \ip s ,a) represent the entangled states described in Section[3j and the splitting 
in their energies is decreased. It can be shown that the polarization of any cell in the unbiased 
line is zero if the line is in the steady state described by Eq. |HJ 

Figures [6] and [7] show the results of two simulations of six cells in a one-dimensional chain, 
with periodic boundary conditions and nearest neighbour coupling only. 7 is modulated between 
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1 meV and 200 meV for the first simulation, and between 1 meV and 1000 meV in the second. 
Ek = 108.5 meV in both simulations. The Hamiltonian in equation [3] is used, so that intercellular 
correlations are completely included. In both simulations, the tunneling barriers in cells 1, 2, and 
3 are initially high (meaning that 7 is low) while the tunneling barriers in cells 4, 5, and 6 are low 
(so that 7 is high). As time progresses, the tunneling barriers in cell 4 are raised, allowing it to 
polarize, and the barriers in cell 1 are lowered, which causes it to depolarize. Next, cell 5 has its 
barriers raised as the barriers in cell 2 are lowered, and so on. Because of the periodic boundary 
conditions, the bit packet moves cyclically through the six cells. Throughout the process, the 
Hamiltonian changes quasi-adiabatically, by which we mean that d^y/dt <C E\jh. For our initial 
state, we create a nearly full negative polarization in the three active cells by taking the normalized 
sum of the two lowest energy eigenstates. The resulting cell polarizations are plotted as a function 
of time with and without dissipation in each case. A third curve plotting the exponential decay, 
proportional to e~*/ r , is also shown in each plot as a reference. 

In the simulation of Figure |6j as time progresses, the bit packet travels along the periodic line, 
and a polarization is maintained when there is no dissipation. Because the Hamiltonian is changed 
quasi-adiabatically, no kinks are created as the bit packet moves, and the bit packet evolves 
qualitatively the same as it would if it were stationary. Because the tunneling barriers are never 
lowered completely (7 has a maximum value of 200 meV), the cells never reach zero polarization. 
The interaction between the three active cells at a given time, plus the residual polarization in the 
"inactive" cells, slows the coherent oscillations to the point where the bit packet is fully polarized 
over the period of the simulation, as long as there is no dissipation. However, the steady state still 
has zero polarization in all cells, so the bit packet loses polarization exponentially when relaxation 
is included in the model. The slow oscillation frequency demonstrated in this simulation simply 
permits the cells to maintain their polarizations over several clock cycles. 

In the simulation of Figure [7j since the tunnelling barriers are lowered enough to completely 
depolarize the cells (7 is now allowed as high as 1000 meV), fewer cells are "on" at any given time, 
and the coherent oscillations from negative to positive polarization state are sufficiently fast to be 
observed over the time scale of the simulation. As in Figure [6j the dissipation of energy does indeed 
bring the system exponentially to its steady state of zero polarization. In this simulation, because 
of the visible coherent oscillations, the polarization does not follow the decaying exponential as it 
does in Figure |6j but instead the amplitude of the oscillation is proportional to the exponential 
decay. 

5 Discussion 

The simulations in Section [4] show that the ICHA does not always provide a good approximation 
to the full quantum mechanical model of an isolated array of QCA cells. In the case of bit packets 
which are not coupled to a driving cell, the exponential relaxation causes any initial polarization 
to be lost. Within the framework of equations [3j [7j and [8j the maximum time for a clocked 
computation will be limited by the loss of classical information either through energy relaxation or 
through oscillations of the type shown in Figures 5 and 7. In principle, coherent oscillations do not 
imply a loss of information (the information can be retrieved by a carefully timed measurement), 
however, the oscillations may be much more rapid than the duration of a measurement, e.g. in 
the case of a small number of cells, and will in general suffer phase decoherence, which will erase 



12 




Figure 6. Simulations of a six cell line with periodic boundary conditions. The value of 7 is adiabat- 
ically raised and then lowered between 1 < 7 < 200 meV, with a period of 217h/Ek, with a different 
phase for each cell. The circled line shows the polarization of each cell in the absence of any energy 
dissipation. The solid line shows the polarization of each cell with an energy relaxation time of 
434// /7T/,.. For reference, a decaying exponential as a function of r is shown with the dashed line. 



even this phase information. All the simulations presented here were done within the two-state 
approximation of equation |3j which is an approximation of equation [Tj As mentioned previously, 
working directly with equation [T] is very computationally expensive. However, based on the 
conceptual understanding outlined in this paper, in addition to other computational results on 
very small groups of cells (not shown here), we believe that these new features of clocked QCA are 
not an artifact of the two-state approximation, but are the result of the inclusion of intercellular 
correlations and the use of the relaxation time approximation. 
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Figure 7. Simulations of a six cell line with periodic boundary conditions. The value of 7 is adiabati- 
cally raised and then lowered between 1 < 7 < 1000 meV, with a period of 217H/E}., with a different 
phase for each cell. The circled line shows the polarization of each cell in the absence of any energy 
dissipation. The solid line shows the polarization of each cell with an energy relaxation time of 
A3Ah/ E k . For reference, a decaying exponential as a function of r is shown with the dashed line. 
Half a period of a coherent oscillation can be seen over the course of the simulation, flipping the 
polarization state from negative to positive. 



Previous simulations have used the ICHA to predict the latching of lines and the successful 
propagation of clocked pulses across QCA arrays (even when tunnelling rates were varied by only 
one order of magnitude pi). These have seen some degree of validation from experimental work 
on lithographically defined QCA-like systems [IO] [ll]. The results presented here show that these 
experiments do not embody the dynamics of equations [3] and [7j The use of intermediate quantum 
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dots and multiple tunnel junctions in these experiments makes it possible to tune the tunneling 
time from Wps to ~ 3000s [10|. This drastic change in tunneling rates amounts to a crossover 
from a regime where quantum tunneling is important to one where electrons are completely 
localized. The normal relaxation dynamics are suppressed because the system is strongly coupled 



to the environment through a z , which has the effect of localizing charge 30 . In molecular and 
atomic implementations, it is likely that the range of tunneling rates that allows this crossover 
from quantum dynamics to classical dynamics by directly modulating tunnel barriers, will not 
be achievable. A more promising approach to achieving such a crossover would be to change the 
internal dynamics of a bit-packet by changing its size. In section |3j we showed that longer lines 
exhibit slower dynamics and greater bistability. This suggests that for sufficiently large bit packets, 
there may be sufficient bistability that the coupling to the environment results in localization of 
charge, even with a limited range of tunneling rates available. The precise tunneling rate and bit 
packet size that will allow localization of charge with sufficiently slow relaxation dynamics will be 
determined by the nature and strength of the coupling of the environment to the QCA system, 
which will in turn compete with the QCA system's internal dynamics. A detailed analysis will 
need to take into account the microscopic details of a specific implementation. 

The simulations performed in this work clearly identify depolarization due to quantum correla- 
tions as a critical issue for classical computation using clocked QCA at the molecular and atomic 
scale. To our knowledge, this phenomenon has not been acknowledged in any of the previous work 
on QCA. However, it is important to note that this loss of polarization occurs only when a set 
of cells becomes isolated from a perturbing influence such as a fixed driver cell, and thus, QCA 
systems operating within a single clocking zone will not experience such an effect. 



6 Conclusion 



In this paper, we have assessed some limitations of conventional approaches to QCA simulation. 
Full quantum mechanical calculations show that the ground state of an unbiased cell, or of a line of 
cells, is a superposition of the two fully aligned states, and thus holds no polarization and carries 
no information. When the assumption of exponential relaxation to a thermal steady state is made, 
we find that cells, or groups of interacting cells, lose their polarization over time unless they are 
influenced by a fixed driver. This is the case even if the cells start with an initial polarization 
before being decoupled from a fixed driver cell. This depolarization effect was not predicted in 
previous QCA simulations, which had often predicted a false latching mechanism among cells that 
would allow them to retain their polarization even in the absence of a driver cell. We have found 
that this failure is related to the ICHA which neglects correlations and shows hysteresis in array 
polarization. Although lithographic QCA systems have managed to avoid this problem thanks 
to their inherently long tunnelling times, the molecular and atomic implementations of QCA 
required for room temperature operation will likely behave in a more purely quantum mechanical 
way, so that the solutions of the many-cell Hamiltonian need to be included. Only an appropriate 
and sufficiently strong interaction of a QCA array with its environment will make clocked QCA 
possible. 

While these findings do not affect the original concept of ground state computing with QCA, 
they may require a reconsideration of QCA architecture, specifically relating to clocking and 
memory devices. Because the simulations presented here are still for a very highly idealized 
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model of QCA behaviour, and ignore, among other things, the specifics of the interaction with 
the environment, we do not claim that clocking or memory are impossible in QCA. Effective 
clocking requires a tuneable change from quantum mechanical behaviour to classical. It remains 
to be shown how and if this can be achieved at the molecular scale. This underscores the need 
for a more sophisticated theory of QCA operation, which should include implementation-specific 
dynamics beyond the phenomenological relaxation time approximation. 
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